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IN A CHARGING ENVIRONMENT 
Carolyn K. Purvis 
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Lewis Research Center 
ABSTRACT 

The sensitivity of predicted equilibrium potential to changes in secondary 
electron yield parameters has been investigated using MATCHG, a simple charging 
code which incorporates the NASCAP material property formulations. It has been 
found that equilibrium potential is a sensitive function of one of the two 
parameters specifying secondary electron yield due to proton impact (®p) and 
of essentially all the parameters specifying yield due to electron impact. It 
is further found that information on the electron generated secondary yield 
parameters can be obtained from monoenergetic beam charging data if charging 
rates as well as equilibrium potentials are accurately recorded. 

INTRODUCTION 

Charging of geosynchronous spacecraft during geomagnetic substorm activity 
is modeled in terms of currents to spacecraft surfaces, with the condition for 
equilibrium being that the net current to a surface element be zero. This net 
current is the sum of incident, emitted and conducted currents, and depends upon 
environment, surface material properties and system electrical and physical con- 
figuration. The environment of concern consists of fluxes of kilovolt electrons 
and ions (H + ) injected into the magnetosphere during substorm activity, plus 
solar photons. Charged particle emission processes generally considered to be 
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important for charging calculations include secondary electron production by 
both ion impact and electron impact, electron backscattering and photoelectron 
emission. System configuration concerns include capacities among various system 
con> r '- or.ents and shadowing of spacecraft surfaces. 

Presently available computer codes have been designed to calculate charg- 
ing, and incorporate algorithms for computing currents due to the important par- 
ticle emission processes. These require as input material property parame- 
ters. Many of the parameters needed to specify electron emission yields are 
poorly known or unknown for common spacecraft surface materials, such as kapton, 
teflon, mylar, solar cell cover slips and paints. An experimental program to 
determine accurately all of the required properties for all materials of inter- 
est would be prohibitive. One approach to reducing the magnitude of the experi- 
mental task is to conduct computer studies to identify the importance of the 
material property input parameters in determining potentials attained by sur- 
faces in kilovolt charged particle environments. Such a study, utilizing the 
emission yield formulations of the NASA Charging Analyzer Program (NASCAP) 

(ref. 1) is reported here. A discussion of how the NASCAP material parameters 
are used to calculate yields is given in reference 1. In what follows, the 
present study is described and parameters varied are identified; results are 
presented; and implications of the results are discussed. 

Two points should be noted at the outset. First, the present study is not 
an evaluation of the material property formulations used in NASCAP; it examines 
the impact of input parameters on predictions based on the NASCAP formulations. 
Second, the study does not predict potentials for any particular material; 
rather it examines the effect on predicted potential of various material prop- 


erty parameters. 
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STUDY DESCRIPTION 

This study focuses on the influence of those parameters which determine 
yields of secondary electrons due to electron and ion impact on the potentials 
attained by surfaces in charging environments. In particular, it examines the 
influence of these parameters on predicted equilibrium potentials of surfaces 
exposed to an isotropic Maxwellian particle distribution and on predicted charg- 
ing behavior of surfaces exposed to a monoenergetic electron be.im. 

The charging calculations have been performed with MATCHG (ref. 2), a code 
which uses the material property formulati'ins of NASCAP to calculate surface 
charging, but does not consider multidimensional effects, photoemission, or 
leakage currents. MATCHG models the response of a "capacitor" with one "plate" 
grounded and the other exposed to either a monodirectional , monoenergetic elec- 
tron beam (essentially, an infinite flat plate approximation) or an isotropic 
Maxwellian flux of electrons and ions (a spherical probe approximation). It is 
therefore a one dimensional analytical model which uses the material property 
parameters to determine net current to the exposed surface. Charging rate in 
this model depends on the specified dielectric constant and thickness of the 
material. Thus, charging calculations done with MATCHG do not give a good re- 
presentation of the charging of a complex spacecraft, but this code does provide 
an efficient means of identifying the influence of material property parameters 
on surface charging. All calculations were performed with yields calculated for 
normally incident primaries (see refs. 1 and 2). 

The Maxwellian environment used was characterized by electron and ion tem- 
peratures of 10 keV and number densities of l/cm^. The electron beam was char- 

2 

acterized by a beam voltage of 10 kV and a current density of 1 nA/cm . 


Calculations of equilibrium potential (Cp ) in the Maxwellian environment 
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and charging history in the beam environment were made for a variety of combina- 
tions of secondary yield parameters. The yield of backscattered electrons, 
which in the NASCAP/MATCHG formulation depends only upon the atomic number ('/,), 
was not varied systematically; however, calculations indicate that the back- 
scatter yield varies slowly with Z. Consequently, cp is not very 
sensitive to variations in Z. For the results presented here, Z was held 
constant at a value of 13. This is the correct value for aluminum, and :'s 
reasonably close to effective Z for polymers of interest. The parameters 
whose effects have been investigated in some detail are those which in the 
NASCAP/MATCHG formulation determine the yields of secondary electrons due to 
electron and ion impact. These are, for ion impact, the yield for 1 keV protons 
incident on the surface (6^) and the primary ion energy at which maximum yield 
is attained (E^). For electron impact, the parameters are the maximum yield 
(6 ), the primary electron energy for maximum yield (E ), two range coefficients 
^ r l ,r 2^ anC * tW ° ex P onents * n 2 ) * Th e four parameters define the elec- 

tron range according to 


n 1 n j 

R = r x E 1 + r 2 E (1) 

where E is the primary electron energy at impact. It is possible to set one 
of the coefficients to zero to obtain a single exponential form for R. This 
has been done for parts of the present study in order to clarify the impacts of 
the four parameters in R. There is also an option in the NASCAP/MATCHG mate- 
rial formulation under which the codes will generate a range coefficient and an 
exponent from the atomic number, atomic weight and material density according to 
formulae due to Feldman (ref. 3). This option has not been used for this study 
because the range exponent in this case is determined from the atomic number 
which, as noted above, also determines backscatter yield. 
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PROTON-GENE KATE1) SECONDARY ELECTRONS 
As noted above, two parameters, 6^ and Ep, determine the yield of 
secondary electrons due to ion impact as a function of ion energy at impact. 
The yield is calculated from 


6(E) 


CE 


1/2 


(1 + E/E ) 
P 


( 2 ) 


where C is a constant numerically equal to 6 for E in keV. These were 

P 

varied systematically, and equilibrium potentials calculated for several values 

of the electron-generated secondary yield parameters 6 and E . The 

m m 

range parameters used are the nominal NASCAP parameters for aluminum. Results 

are illustrated in figure 1. The figure shows predicted equilibrium potential 

(cp ) as a function of 6 parameterized by E and 6 . The curves 
v eq' p p m 

labeled "6^ = 0" indicate potentials calculated with secondary electron emis- 
sion due to electron impact set to zero. Thus they represent an upper limit 
on cp for the given Z and environment, independent of 6^, E m and 
the range parameters. It is clear from these curves that cp £ s very sensi- 
tive to changes in 6^, particularly when 6^ is small, but relatively insen- 
sitive to changes in E p . This result holds true for the more realistic 
cases, 6 m t 0, as is illustrated in the figure. It is expected that the basic 

result that cp is sensitive to 6 and insensitive to E„ also holds 
eq p p 

for other Maxwellian environments so long as kT^ < E^. if the ion tempera- 
ture were to exceed E cp would be more sensitive to E_. However, 
since E^ is usually expected to be about 40 to 100 keV, and ion temperatures 
at geosynchronous are generally expected to be less than 40 keV, it seems rea- 


sonable to expect that F is not a critical parameter for charging calcula- 
tions. b on the other hand is a critical parameter. 

P 
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ELECTRON GENERATED SECONDARY ELECTRONS 
Yield of secondary electrons due to electron impact is determined by a 
total of six parameters: 6^, the maximum yield; E m , the primary energy for 

maximum yield; and the four range parameters described above. The yield is cal 
culated from 


fi e (E) 



dx 


(3) 


Here, |dE/dx! is the energy loss rate of the primary electron in the material, 
dx is an element of path length, 8 is the angle of incidence of the primary 
electron, R u is the stopping distance, and C^ and C£ are constants. 

The energy loss rate is related to the range R by 


dE 

dx 




(A) 


with Ej the energy of the incident primary at impact. NASCAP takes the 
first two terms of a Taylor series expansion, and uses 



(5) 


with the empirical formula 

n l n 2 

R = r x E 1 + r 2 E ( 1) 

as noted above. 

The stopping distance is defined from 

/'" ia 

0 


dx = E 


( 6 ) 
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With these expressions for R u and |dE/dx|, the constants Cj^ and C 2 

in equation (3) are determined by requiring that G g (E) have the user- 

specified values of 6 and E . 

m tn 

Potentials in Maxwellian Environment 

Results of varying 6^ and E ]n for a fixed set of range parameters 

are illustrated in figure 2. This figure shows cp as a function of 6 

eq m 

parameterized by E and 6 in a format analagous to figure 1 for the 

Fn p 

proton yield parameters. The range parameters are the same as those used for 

the calculations summarized in figure 1. The "6 = 0" curves again represent 

P 

an artifical upper limit to q> for a given 6 ,E„ combination and the 

eq b m* m 

stated values for 2 and range. The curves in figure 2 indicate that <p is 
quite sensitive to 6 , and somewhat less sensitive to E„. The sensitivity 
of $ to both these parameters is more pronounced for values of 6^ > 1. 

The discussion to this point has, as noted above, considered a particular 
combination of range parameters. It would be convenient if these parameters did 
not have a significant effect on <p . Unfortunately, their effect is pro- 
found. In order to understand the effect of the various parameters on , 

it is useful to examine their effect on the secondary yield curve itself. For 
the purpose of examining effects on charging behavior, it is helpful to plot 
"normalized" yield curves, such as the one shown in figure 3. This curve shows 
the yield of secondary electrons due to electron impact normalized to 6 

m 

(i.e., 6/6 ) plotted as a function of electron energy at impact normalized to 

E (i.e., E/E ). The shape of this normalized curve depends only upon the 
m m 

range parameters. The roles of the various range parameters in determining the 
shape of this normalized curve can be clarified by considering first a single 
exponential expression for the range (i.e., setting one of the coefficients to 
zero) . Then 
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rE 


(7) 


In this case, varying the value of r was found to have no impact on the yield 
curve; varying n had a profound effect. Figure 4 shows normalized yield cur** 
ves for four different values of n for E/E m greater than 1, which is the 
ranje of interest for spacecraft charging calculations. The yield is clearly 
very sensitive to n. It is thus expected that 9 is very sensitive to n. 

The yield curve calculated using a double exponential range expression, 

R » rjE n l + r^ n 2, must be between the two "single exponential range" cur- 
ves determined by n^ and n 2 < The exact shape of the yield curve depends 
on both exponents and on the ratio of '.he coefficients, r^/^; i*- does not 
depend on the magnitudes of r^ and Figure 5 illustrates the effect 

of the coefficient ratio on predicted yield for the case R ■ r^E^*^ + 

The single exponent curves for R ■ rE^'^ and R * rE^*® are 
drawn in for comparison. The yield curve falls off more slowly with increasing 
energy as r \l^2 increases. The differences in yield between the double 
exponential ranges and the single exponential range with n = 2 are smaller for 
large energies. The sensitivity of to r^/r 2 is thus expected to 

depend upon the temperature of the environment relative to E m » For example, 
for the 10 keV Maxwellian environment being used here, and with E m =0.3 keV, 
kT e /E m ■ 33 1/3. The curves in figure 5 indicate that, at E/E m = 33 1/3, 
the yield for x \l^2 = 0.25 is very similar to that for n M 2; the yield 
for r^/r 9 = 1 is somewhat larger; and that for = ** i® signifi- 

cant Iv larger. It is thus expected that cp will be reduced in magnitude 

eq 

a s r ^ / r , is increased. 
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The dependence of 9 on n for Che single exponential range expres- 

eq 

•ion, and on n^ (the lower exponent) and r^/^ f° r values of n 2 of 

1.6 and 2.0 is shown in figure 6. The other secondary yield parameters 

(6 ,E_, etc.) were held fixed at the values indicated on the figure. It is 
b m* 

clear from this figure that <p is very sensitive to n in the single expo** 
nential range caue. In the double exponent case, 9 depends on all three 
variables, i.e., on nj and rj/r^, the sensitivity of <p to 
these variables depends on the values of the variables. 9 is very sensi- 
tive to nj’, its sensitivity to n^ depends on f]/ r 2 > being high when 
r^/r 2 is large. 

Comparison of figures 4, 5, and 6 indicates that for a particular combina- 
tion of the other parameters and environment, it is possible to estimate Cp^ 
for double exponential range expressions from values of <p calculated from 
single range expressions by comparing the yields at kT e /E m . This provides 
a rough estimate rather than a precise one because the calculation of current 
due to secondary electronics from electronic impact involves integration of the 
secondary yield curve multiplied by the distribution function f(E) of incoming 
electrons, i.e., 

* Const y* 6 e (E)f(E) dE (8) 


Thus this current, which is an important component in t lie current balance deter- 


mining Cp , depends on the values of the yield over a range of energies, not 
eq 

simply that at kT e /E m . 

The dependence of Cp on 6 , E and n is summarized in figure 7 , 

eq m m 

where cp is plotted as a function of 6 , parameterized by n and E . 

sq m m 

The curves in figure 7 indicate that both the range exponent and 6 m are 

critical parameter in determining CP , with E of lesser importance. 
r eq m 
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The range parameters taken collectively, along with 6 , are the critical 

m 

determiners of ip for double exponential range expressions. The im- 
portance of the individual parameters In determining *|> depends on their 
values, as discussed above. 

Charging in Electron Bern 

Spacecraft are charged by distributions of particles, which are usually 
modeled as Maxwellian distributions. In contrast, laboratory studies of charg- 
ing behavior of materials have generally been performed using monoenerget ic 
electron beams as particle sources. It is therefore of interest to understand 
how charging of material samples in a monoenerget ic beam depends upon the mate- 
rial property parameters in the models and how laboratory charging data may be 
used to infer parameter values where direct measurements of secondary yield 
against energy are not available. 

Figure 8 shows MATCHC results for <p as a function of 8 parame- 
w 6() in 

terized by E m and n using a 10 kV beam of electrons as the environment. 

Again, cp depends on all of these variables. Comparison with figure 7 in- 
eq 

dicates that <P in the beam environment has a different dependence on n 
and 6^ than in a Maxwellian environment, particularly for small In 


fact, cp in a monoenergetic beam environment is determined by 

eq 

,<V 8 - 1*^1) ■ E n 


(9) 


where is the beam voltage, e the magnitude of the electronic charge and 

Ejj the second unity crossover of the total yield curve (i.e., true secondary 
yield 6^ plus backscatter yield). This is exactly true in the MATCHG ap- 
proximation of an infinite flat plate geometry and no leakage, and approximately 
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true in tile three dimensional NASCAP calculation* and in the true laboratory 
situation for high resistivity samples. 

This implies that measurements of equilibrium potential in electron beams 
essentially give information about only one point on the yield curve. While 
Ejj is an important point on the yield curve, knowing it is insufficient. As 
has been noted, equilibrium potential in a Maxwellian environment depends on the 
shape of the entire yield curve, multiplied by the environmental electron dis- 
tribution function and integrated over energy (eq. (8)). Even if backscatter 
yield and 6 ffi and E n are known, E^^ does not uniquely determine the 
range parameters except in the case of a single exponential range expression. 

While in the monoenergetic beam case depends only on Ejj, the 

charging behavior as a function of time depends upon the shape of the total 
yield curve at impact energies greater than E n . Charging rate depends upon 
capacitance and charging current. The charging current depends on the beam cur- 
rent and the emitted current. The emitted current, in the MATCHG approximation, 
is determined by the value of the total yield at an impact energy 

E. - e(V g - |*P| ) (10) 

where |^l is the magnitude of the surface potential. Large yields result in 
small net currents and therefore in low charging rates. 

Figure 9 shows predicted charging behavior for 6^ * 3, g m ■ o.3 keV 
and several range expressions for a 10 kV beam. Referring back to figures 4 and 
5, one can see that the differences in charging behavior in figure 9 reflect the 
shapes of the secondary yield curves. 

The initial charging rate for a sample with fixed capacitance is determined 
by the total yield at the be«an energy eV g . It is thus expected to depend on 
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the various secondary yield parameters, and on beam energy. The dependence of 
initial charging rate in a 1U kV beam on the secondary yield parameters is il- 
lustrated in figure 10. Thia figure indicates that initial charging rate is a 
very sensitive function of n for all values of 6^. 

Thus, accurate data on charging behavior versus time can be used to iden- 
tify suitable values for the range parameters if 6 ffi * n d E m are known. 

SUMMARY AND CONCLUSIONS 


Using the NASCAP /MATCHC material property formulations, the predicted 
equilibrium potential of a surface in an isotropic Maxwellian plasma has been 
found to be sensitive to almost all of the secondary electron yield parameters. 

The yield of secondary electrons due to proton impact is dependent on two 
parameters, the yield for 1 keV incident protons (6 p ) an d the energy for maxi- 
mum yield (E ). 9. has been found to be very sensitive to 6 , parti- 
p eq ' p r 

cularly for values of 6 less than unity. Accurate values of 6 are 

P P 

thus required if surface potential* in s^ace are to be predicted accurately. 

This is of particular concern because in general 6 has not been measured 

P 

for spacecraft aurface materials. In addition, laboratory investigations of 
charging reported to date have used electron beams and thus do not provide any 
information on ion-generated secondary electrons. 9 has been found to be 
rather insensitive to the energy parameter E p ; it is thus not critical to 
know the value of this parameter accurately, provided that K greater than 
the temperature of the ion distribution. 

Secondary electron yield due to electron impact is determined by a total of 
six parameters: 6^, E m , rp r 2 , n^ and n 2 , where the last four 

are coeffieie*** and exponents in the range expression given by equation (1). 

One of the coefficients can bo set equal to zero, resulting in a range expres- 
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sion of the form R ** rE n . The sensitivity of <p to the various parame- 

ters in an isotropic Maxwellian environment has been investigated using both 
types of range expression. 

Using single exponential range expressions, it was found that cp^ was 

very sensitive to both n and 6 , somewhat less sensitive to E„, and com- 

pletely insensitive to r. The degree of sensitivity to small changes in any 

one of the parameters (except r) was found to depend upon the values of all of 

them. In general, n and 6 ffl were critical parameters, with 6^ becoming 

more critical as 6 was increased or n decreased. The value of E_ 
m m 

becomes more important in determining cp as n and 6 increase. 

eq m 

Calculations using double exponential range expressions indicated that 

" was sensitive to n, , n 0 , r,/r,. 6 , and E , but not to the 

eq 1’ 2’ 1' 2* m m* 

magnitudes of r^ and ^ individually. Again, the degree of sensitivity 
of 9 to variations in any one of these parameters depended on the values 
of them all. For small values of 2 , (i.e., < 1), the important parame- 
ters are 6^, and E m , as in the single exponential range case. For 

larger values of r^/r £, n^ becomes increasingly important in deter- 
mining cp 

eq 

Thus it is necessary to obtain reasonably accurate values for all the 

parameters characterizing secondary yield due to electron impact, except for the 

magnitudes of the range coefficients r^ and ^ For some materials of 

interest to spacecraft charging such as teflon, kapton, and mylar (ref. 4), data 

on secondary yield as a function of impact energy, and thus values for 6^ 

and E , are available. For these materials it is possible to infer the 
m 

needed range parameters from careful charging history data, since both the 
equilibrium potentials and the charging rates in monoenergetic electron 


beams are sensitive to the secondary yield parameters. Charging rate data 
is essential to identification of suitable parameter values, because it 
can be used to infer yield as a function of primary energy at impact. 

Data on equilibrium potentials in monoenergetic beams essentially identi- 
fies only the energy at which the total yield is equal to one. 
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